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ABSTRACT 

Accelerating expansion of the Universe is a great challenge for both physics and 
cosmology. In light of lacking the convincing theoretical explanation, an effective de- 
scription of this phenomenon in terms of cosmic equation of state turns out useful. 

The strength of modern cosmology lies in consistency across independent, often 
unrelated pieces of evidence. Therefore, every alternative method of restricting cosmic 
equation of state is important. Strongly gravitationally lensed quasar-galaxy systems 
create such new opportunity by combining stellar kinematics (central velocity dis- 
persion measurements) with lensing geometry (Einstein radius determination form 
position of images). 

In this paper we apply such method to a combined data sets from SLACS and 
LSD surveys of gravitational lenses. In result we obtain the cosmic equation of state 
parameters, which generally agree with results already known in the literature. This 
demonstrates that the method can be further used on larger samples obtained in 
the future. Independently noticed systematic deviation between fits done on standard 
candles and standard rulers is revealed in our findings. We also identify an impor- 
tant selection effect crucial to our method associated with geometric configuration 
of the lensing system along line of sight, which may have consequences for sample 
construction from the future lensing surveys. 
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1 INTRODUCTION 

The present acceleration of the cosmic expansion is a funda- 
mental challenge to standard models of both particle physics 
and cosmology. Discovery of this phenomenon on the Hubble 
diagrams obtained from the SNIa surveys (Riess et al. 1998, 
Perlmutter et al. 1999, Riess et al. 2004, Wood- Vasey et al. 
2007; Davis et al. 2007; Kowalski et al.2008 ) in combination 
with independent estimates of the amount of baryons and 
cold dark matter (Spergel et al. 2003, Eisenstein et al. 2005) 
led us to believe that most of the energy in the Universe 
exists in the form of mysterious dark energy. 

The new physics of dark energy may lie in the na- 
ture of gravity, the quantum vacuum, or extra dimensions. 
Concerning the first possibility there exists an increasing 
body of literature (e.g. Buchert 2001, Rasanen 2004, Ellis & 
Buchert 2005, Wiltshire 2007) pointing out that if one at- 
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tempts to average out local sources of gravity (galaxies and 
clusters) in order to obtain the smoothed description of the 
Universe in the large such averaging procedure not com- 
muting with temporal evolution could manifest as an addi- 
tional source term in the energy-momentum tensor. Within 
the second possibility our ideas about the quantum vacuum 
are expressed by either introducing cosmological constant A 
or some time evolving scalar field (quintessence). The last 
possibility is to contemplate modifications to the Friedman- 
Robertson- Walker models arising in brane-world scenarios. 
Irrespective of the theoretical approach chosen a common 
point with the observations usually occurs at the level of the 
w{z) coefficient in an effective equation of state p = w{z)p 
for dark energy. 

The potential of constraining dark energy models with 
SNIa data, even though ever increasing, would not be suffi- 
cient if taken alone in separation form the other approaches. 
Indeed, the power of modern cosmology lies in building up 
consistency rather than in single, precise, crucial experi- 
ments. Therefore, every alternative method of restricting 
cosmological parameters is desired. In this spirit a number 
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of combined analyses involving Icnsing statistics (Silva & 
Bcrtolami 2003), CMBR measurements (Spergel et al. 2003, 
Hinshaw et al. 2009), age-redshift relation (Alcaniz, Jain & 
Dev 2003), x-ray luminosities of galaxy clusters (Allen et al. 
2008) or the largo scale structure considerations (Eiscnstcin 
et al. 2005) have been performed in the literature (references 
above being far from complete). 

In this paper we use strongly gravitationally lensed sys- 
tems for providing additional constraints on dark energy 
models. The idea of using such systems for measuring the 
cosmic equation of state weis discussed in Biesiada (2006) 
and also in more recent paper by Grille et al. (2008). The 
first (to our knowledge) formulations of this approach can 
be traced back to Futamasc & Yoshida (2001). Next sec- 
tions outline the method, the sample used and cosmological 
scenarios tested. The last section presents the results and 
conclusions. 



2 THE METHOD 

Strong gravitational lensing occurs whenever the source, the 
lens and observer are so well aligned that the observer- 
source direction lies inside the so called Einstein ring of the 
lens. In cosmological context the source is usually a quasar 
with a galaxy acting as the lens. Although strong lensing 
by clusters is known and number of such cases increases, we 
will be concerned with galaxies acting as lenses. For detailed 
theory of gravitational lensing see e.g. Schneider, Ehlers & 
Falco (1992). Strong lensing reveals itself as multiple im- 
ages of the source. The image separations in the system de- 
pend on angular-diameter distances to the lens and to the 
source, which in turn are determined by background cosmol- 
ogy. This opens a possibility to constraining the cosmologi- 
cal model provided that we have good knowledge of the lens 
model. 

Since the discovery of the first gravitational lens the 
number of strongly lensed systems increased to a hundred 

(in the CASTLES database §) and is steadily increasing fol- 
lowing the new surveys like SLAGS (Sloan Lens AGS Survey 

^). It turns out that in vast majority of cases the lens is a 
late type E/SO galaxy. This could be understood since ellip- 
ticals being a latecomers in hierarchical structure formation 

arc created in mergers of low-mass spiral galaxies. Hence 
they are more massive than spirals and because the Ein- 
stein ring radius scales with mass, the probability of their 
acting as lenses is higher. Despite they lack bright kinematic 
tracers at large radii (e.g. like HI in disk galaxies) and thus 
their kinematics is more difficult to measure, there exists an 
increasing evidence that their mass density profile can well 
be approximated by singular isothermal sphere (SIS) model 
(or a variant thereof called singular isothermal ellipsoid - 
SIE). 

Now, the idea is that formula for the Einstein radius in 
a SIS lens (or its SIE equivalent) 

e^_4n^— (1) 

§ http: // www. cfa. harvard . edu/castles / 
^ http://www.slacs.org/ 



depends on the cosmological model through the ratio of 
(angular-diameter) distances between lens and source and 
between observer and lens. The angular diameter distance 
in flat Friedmann-Robertson- Walker cosmology reads: 

N 1 c dz' ,„, 

where Ho is the present value of the Hubble function and 
h{z; p) is a dimensionless expansion rate dependent on red- 
shift z and cosmological model p arameters p. For example 
in flat AGDM model h{z; p) = ■\/f2m,(l + z)^ + SIa we have 
Qa = 1 — hence p = {f^m} in this case. Expansion 
rates in other cosmological scenarios are given in Section 4. 
From now on we will assume spatial flatness of the Universe 
since it is strongly supported by independent and precise ex- 
periments, e.g. a combined WMAP5, BAG and supernova 
data analysis gives Qtot = l.OOSOlg QQg? (Hinshaw et al. 
2009). The sample upon which wc work is small, and ad- 
dition of (otherwise well constrained) curvature parameter 
would only distort the results. 

Provided one has reliable knowledge about the lensing 
system: i.e. the Einstein radius 6e (from image astrometry) 
and stellar velocity dispersion asis (form central velocity 
dispersion ao obtained from spectroscopy) one can use it 
to test the background cosmology. This method is indepen- 
dent on the Hubble constant value (which gets cancelled in 
the distance ratio) and is not affected by dust absorption 
or source evolutionary effects. It depends, however, on the 
reliability of lens modelling (e.g. SIS or SIE assumption) 
and measurements of ao- Hopefully, starting with the Lens 
Structure and Dynamics (LSD) survey and the more recent 
SLAGS survey spectroscopic data for central parts of lens 
galaxies became available allowing to assess their central ve- 
locity dispersions. In practice central velocity dispersion ao 
is estimated from the velocity dispersion within Re/ 8 where 
Re is optical effective radius. Thorough discussion of these 
issues can be found in (Treu et al. 2006, Grillo st al. 2008) 
where the arguments in favor of using ao as representative 
to asis are presented. Moreover, there is a growing evidence 
for homologous structure of late type galaxies (Treu et al. 
2006, Koopmans et al.2006, 2009) supporting reliability of 
SIS/SIE assumption. In particular it was shown in (Koop- 
mans et al. 2009) that inside one effective radius massive 
elliptical galaxies are kinematically indistinguishable from 
an isothermal ellipsoid. 

In the method used in this paper cosmological model 
enters not through a distance measure directly, but rather 
through a distance ratio 

^ - ^ J . - .z. g , l^j 

Jzi hiz';p) 

and respective observable counterpart reads: 
^obs _ 47ra-o 

This has certain consequences both advantageous and disad- 
vantageous. The positive side is that the Hubble constant Ho 
gets cancelled, hence it does not introduce any uncertainty 
to the results. Ou the other hand we have a disadvantage 
that the power of estimating is poor (which could be 
seen by inspection into specific formulae for h{z; p) - see Ta^ 
ble 1 below). Therefore we only attempted to fit ilm in the 
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Table 1. Combined SLAGS + LSD lens sample. The restricted 
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case of ACDM model (where it is the only free parameter in 
flat cosmology) and it was successful only for the restricted 
sample (see below). In other cases we assumed fixed values 
for fim. Cosmological model parameters (coefficients in the 
equation of state) have been estimated by minimizing the 
chi-square: 

(4) 

where the sum is over the sample and crl) j denotes the vari- 
ance of X>°^^ (contextual use of the same symbol for variances 
and velocity dispersions should not lead to confusion). In cal- 
culating gt> we assumed that only velocity dispersion errors 
contribute and the Einstein radii are determined accurately. 



3 SAMPLES USED 

We used a combined sample of n = 20 strong lensing sys- 
tems with good spectroscopic measurements of central dis- 
persions from the SLACS and LSD surveys (essentially the 
same sample as used by Grillo et al. (2008)). Original data 
concerning SLACS sample came from Treu et al. (2006) (see 
also an erratum (Treu et al. 2006a) - very important one). 
Data concerning LSD lenses are taken after Treu and Koop- 
mans (2004), Koopmans and Treu (2003, 2002). 

As already noticed in Treu et al. (2006) the SLACS sam- 
ple has an average Dis/Ds ratio equal to 0.58 with an rms 
scatter 0.15. Whereas for their purpose it was advantageous, 
in our context it weakens the performance of the method. 
Therefore we selected a sub-sample of n = 7 lenses with the 
TD ratio deviating from the mean more than rms in either 
direction. It is summarized in Table 1 where the names of 
lenses in the restricted sample arc given in bold. 

For comparison of our results with the data which trig- 
gered the dark energy problem, we also performed fits to the 
SNIa data (n=307 supernovae) using Union08 compilation 



by Kowalski et al. (2008). The Q,m = 0.27 prior was used 
throughout, except in the ACDM model where the fit was 
attempted. 



4 COSMOLOGICAL MODELS TESTED 

Several scenarios have been put forward as an explanation of 
presently accelerating expansion of the Universe. The most 
obvious candidate is the cosmological constant A represent- 
ing the energy of the vacuum. Corresponding cosmological 
model, which turned out to be in agreement with all exist- 
ing (independent and alternative) observations is the ACDM 
model. It is equivalent to w = — 1 in the cosmic equation 
of state p = wp and the only free parameter here is the 
f2m representing the density of baryonic plus cold dark mat- 
ter as a fraction of critical density (as already said spatial 
flatness is assumed). On one hand it is therefore the most 
parsimonious one, but well known fine tuning problems led 
many people to seek beyond the A framework and to develop 
the concept of quintessence. Usually the quintessence is de- 
scribed in a phenomenological manner, as a scalar field with 
an appropriate potential. In first approximation it could be 
tested observationally by promoting w to the role of a free 
parameter to be fitted from the data. However there is no 
a priori reason to expect that w should then be a constant. 
The parametrization of w{z) = wo + 'Waj^ developed by 
Chevalier & Polarski (2001) and Linder (2003) turned out 
to be well suited and robust for such case. In the past, alter- 
native parametrization w{z) = wo + wiz was used (which is 
a truncated Taylor series representation of 'w{z)). Chevalier- 
Polarski-Linder parametrization instead uses an expansion 
with respect to the physical degree of freedom i.e. the scale 
factor (expanded around its present value). Dimensionless 
(i.e. with Ho factored out) expansion rates for respective 
models are given in Table 2. 

For comparison we also performed fits of the models 
considered above to the supernova la data with the same 
prior assumptions (spatial flatness of the Universe and fim,)- 
We have taken Union08 compilation (Kowalski et al. 2008) 
and instead of straightforward fitting rn{z) vs. Dl{z) 
(where m is visible magnitude Dl denotes luminosity dis- 
tance) we used a well known modified approach equivalent to 
marginalizing over the intercept (Nesseris & Perivolaropou- 
los 2005). 



5 RESULTS AND CONCLUSIONS 

Performing fits of different cosmological scenarios (shown 
in Table 2) on the fuU SLACS + LSD sample of n = 20 
strong lensing systems we obtained the equation of state 
parameters displayed in Table 3. In ACDM model, where 
Qrn was the only free parameter we were not able to make a 
reliable fit on the sample considered. As already mentioned 
the reason of this is twofold. First, theoretical observable 
was the ratio of two integrals differing only by the limits 
of integration. Second, in the full sample the Dis/Ds ratio 
is concentrated around a central value of 0.54. Therefore, 
in case of 'simple' dependence (just a factor) of h{z) on a 
parameter (as is the case for Qrn) the bulk of the sample 
only introduces an undesired scatter. More 'sophisticated' 
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Table 2. Equation of state and expansion rates h(z) in the models tested. 

Model EOS p = wp Cosmological expansion rate h{z) 

ACDM w= -1 h{z) = sJO.m (1 + + Ha 

Quintessence w = const. h{z) = ^/^m (1 + z)^ + (1 + 2)3(1+'") 

Chevalier-Polarski-Linder w{z) =wo + Wajf^ h{z) = ^Clm (1 + z)^ + Qq (1 + z)Hi+'^o+wa) exp{^^^^) 



Figure 1. Best fits (dots) and (68 %, 95%) confidence regions for 
Chevalicr-Linder-Polarski parameters in cosmic equation of state 
obtained from the full SLACS+LSD sample of lenses and UnionOS 
SNIa data. A dot corresponding to ACDM model is added for 
reference. 



Figure 2. Best fits (dots) and (68 %, 95%) confidence regions for 
Chevalicr-Linder-Polarski parameters in cosmic equation of state 
obtained from restricted sample of lenses and UnionOS SNIa data. 
A dot corresponding to ACDM model is added for reference. 



dependence (exponent of the integration variable) on the 
equation of state parameters makes it possible to obtain fits 
on this sample. 

For comparison we also report (Table 4) values of these 
parameters best fitted to the UnionOS SNIa compilation 
(Kowalski et al. 2008). One can see that w coefficient ob- 
tained from the full strong lensing sample agrees with re- 
spective value derived from supernovae data (almost whole 
2a confidence interval for w from UnionOS data set lies 



within 1(7 CI from lenses). The value inferred is also in agree- 
ment with the WMAP5 results presented in (Hinshaw et 
al. 2009) including also combined WMAP5, BAG and SNIa 
analysis. Note that this is also consistent with the ACDM 
model. Concerning evolving equation of state in Chevalier- 
Polarski-Linder parametrization confidence regions in the 
{wo,Wa) plane arc shown in Figure 1. It can be seen that 
the concordance model (ACDM) while consistent with SNIa 
data (at 2a level) is inconsistent with strong lensing systems 
method applied here. SNIa results and strong lensing results 
are marginally consistent with each other. Most probably 
this is due to small sample of strong lenses combined with 
systematics discussed above. 

Working on the restricted sample (containing lenses 

with Dis/Ds ratio greater than rms spread around mean 
value) although decreased sample size dramatically (down 
to n = 7) allowed to obtain fits on in ACDM (Table 
3) which turned out to agree with SNIa fits (Table 4) and 
WMAP5 data (Hinshaw et al. 2009). Although the best fit 
for w parameter quintessence scenario is higher than inferred 
from SNIa or WMAP5, yet the 2a interval for UnionOS data 
falls within 2a interval from lenses. Hence the agreement is 
quite good. Similarly fits for wo and Wa are improved (even 
though confidence regions get larger). 

One should note however that systematic shift down- 
wards in {wo, Wa) plane persists. Such shift in best fit- 
ted parameters inferred from supernovae (standard can- 
dles, sensitive to luminosity distance) and BAO or acoustic 
peaks (standard rulers, sensitive to angular diameter dis- 
tance) has already been noticed and discussed by Lazkoz, 
Nesseris & Perivolaropoulos (2007) and by Linder & Roberts 
(2008). Bearing in mind similar mutual inconsistency in the 
Hubble constant values inferred from lensing time delays 
Ho = 52 ±6 km s'^ Mpc~^ (Kochanek & Schechter 2003) 
and from the HST Key Project ffo = 72 ± S fem s"^ Mpc~^ 
(Freedman et al. 2001) our result suggests the need for tak- 
ing a closer look at compatibility of results derived by using 
angular diameter distances and luminosity distances respec- 
tively. It is also worth noticing that the ideas of testing the 
Etherington reciprocity relation between these two distance 
measures have been discussed by Basset & Kuntz (2004) and 
by Uzan, Aghanim & MelUer (2005). 

In conclusion our results demonstrated that the method 
discussed in Biesiada (2006) and extensively investigated by 
Grillo et al. (2008) on simulated data can be used in prac- 
tice to constrain cosmological models. It turned out to give 
reasonable results on already accessible samples of strongly 
lensed systems. Besides the uncertainties related to velocity 
dispersion measurements and their conversion to relevant 
lens model parameters (as well as the impact of SIS as- 
sumption) the issue of systematics associated with Dis/Ds 
ratio in the sample turned out to be important. In partic- 
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ular it implies that strong lensing survey strategies like the 
one adopted in SLAGS survey are better from this point of 
view. Lensing systems are gathered around something like 
0.58 in distance ratio because it is roughly the configura- 
tion for which lensing probability (for a given lens mass) is 
the highest. Earlier searches were focused on source popula- 
tion (quasars) seeking for close pairs or multiples and check- 
ing if they are multiple images of a single source lensed by 
an intervening galaxy. Therefore a high lensing probability 
was an important selection factor there. On the other hand 
SLACS survey is focused on possible lens population (mas- 
sive elipticals) with good spectroscopic data. Using SDSS 
templates spectra arc carefully checked for residual emission 
(at least three distinct common atomic transitions) coming 
from higher rcd.sliifts. Such candidates undergo image pro- 
cessing by subtracting parametrized brightness distribution 
typical for early type galaxies in order to reveal multiple 
images of the quasar. Details can be found in Bolton et 
al. (2006). Therefore, besides the obvious bonus of having 
central velocity dispersion measured, such strategy is bet- 
ter suited for discovering systems with larger Dis/Da ratios 
which in turn can be used for testing cosmological models. 

Finally, one important effect - neglected here - should 
be mentioned, which is the influence of line of sight mass 
contamination. The debate on this issue started with Bar- 
Kana (1996) and Keeton et al. (1997) who were among the 
first to convincingly demonstrate that the effect of large scale 
structure on strong lensing could be significant. More recent 
results on this issue can be found in Dalai et al. (2005) (in the 
context of cluster lensing) or Momcheva et al.(2006). This 
rises the issue of an impact this effect might have on our 
results, since the sample was small. Straightforward naive 
first guess (based on Poissonian statistics) might suggest 
that sample size of order of a few hundred lenses might 
reduce line of sight 'noise' contamination down to a few 
percent. This is however not that simple since the line of 
sight contamination is in fact a systematic effect. Namely, 
massive eaxly type galaxies (ie. typical lenses) prefer over- 
dense environments, so one consistent approach would be 
to follow light rays (ray-shooting simulation) through many 
lens planes (obtained from cosmological N-body simulation) 
up to high source redshift. This was done by Wambsganss 
et al. (2005) with the result that up to = 1 most (i.e. 
95% ) of lenses involved only a single mass concentration, 
whereas for sources at Zs = 3.8 important contribution of in- 
tervening mass could be significant in 38% of strong lensing 
systems. This result suggests that the line of sight contam- 
ination should be addressed separately for each particular 
survey. For the SLACS survey (where the bulk of our sample 
comes from) this was assessed in Treu et al. (2009) where the 
authors found that SLACS lens galaxies are unbiased popu- 
lation (i.e. with environmental effects typical to the over-all 
population of early type galaxies) and typical contribution 
from external mass distribution is small - no more than a 
few percent. Fortunately, the SLACS survey is ongoing i.e. 
the sample of spectroscopically investigated strong lenses is 
growing. However this survey is relatively shallow, so for 
cosmological applications one is forced to combine it with 
deeper surveys (with different designs - hence different sys- 
tematics) and the problem of line of sight contamination 
remains challenging. 
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Table 3. Fits to different cosmological models from SLACS -|- 
LSD lens samples (n = 15 -t- 5 = 20). Fixed value of Clm = 0.27 
was assumed. 



Cosmological model Best fit parameters (with la) X^/'^of 

ACDM not possible 

Quintessence w = -0.9829 ± 0.2415 3.41 

Chevalier-Linder-Polarski = 1.2605 ± 0.8177 3.05 

Wa = -9.4443 ± 4.4193 



Table 4. Fits to different cosmological models from restricted 
SLACS + LSD lens sample (= 7). Fixed value of ftm. = 0.27 
assumed, besides ACDM , where the fit was successful. 



Cosmological model Best fit parameters (with la) /dof 

ACDM Clrn = 0.2660 it 0.2796 1.76 

Quintessence w = -0.6320 ± 0.4461 3.91 

Chevalier-Linder-Polarski wq = 0.3588 ± 1.2453 1.88 

Wa = -3.6301 ± 5.3278 
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